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SPECTRAL  PROPERTIES  OF  PERIODICALLY 
CORRELATED  RANDOM  PROCESSES 


Dr agan ,  Ya . P . 

Physicomechanical  Institute,  AS  UkSSR 
Submitted  20  Jun  70 


In  connection  with  the  importance  of  the  class  of  periodically 
correlated  ( PK )  random  processes  (SP)  in  the  investigation  of  rhythmic 
phenomena  [1]  and  the  insufficient  degree  of  development  of  their 
theory,  a  need  arises  for  a  comprehensive  study  of  them.  In  this 
article,  developing  the  idea  in  work  (2],  the  spectral  properties  of 
PKSP  are  investigated.  By  periodic  correlation  we  understand  the 
periodic  nonstat ionar i ty  of  the  second  order. 

As  it  follows  from  the  results  of  work  [2]  and  the  definition  of 
current  spectrum  [5,  6],  the  function  of  covariation  of  PKSP 

is  presented  in  the  form 

(1) 


in  this  case  the  dependence  of  spectral  function  on  time  is  periodic, 
and  this  function  possesses  the  expansion 


OB  ...  . 


(2) 


where  T  -  period  of  correlation  of  the  process,  and  the  functions 
are  determined  by  the  correlations 


(3) 


1 


pe £ i n i t i on  1 .  As  the  precise  period  of  a  periodic  function  f(') 
we  shall  call  that  least  non-negative  number  T,  for  which  at  all  t  the 
equality  /{t+T)=  /(t)  is  valid.  The  constant  variable  is  a  degenerate 
periodic  function,  since  it  is  periodic  in  the  case  of  any  T  and  in  the 
meaning  of  definition  1  it  does  not  possess  a  precise  period  which  is 
different  from  zero. 

Definition  2.  We  will  call  the  functions  the  correlation 

components  (KK). 

Taking  into  account  that  the  dispersion  of  the  values  of  the 
process 


o*(0  =  6{/.  0). 


(4) 


on  the  basis  of  what  was  said  it  is  easy  to  be  certain  of  the  validity 
of  the  assertion  *. 

*  this  article  for  simplicity  wt  will  asiuiiuc-  lli.'t  i  h<  iii.  lln  i.  *  li  .  i 

expectation  of  the  process  is  equal  to  ^oco.  [end  of  footnote.] 

Theorem  1.  The  current  spectral  function  and  the  dispersion  of 
the  values  of  PKSP  are  periodic  functions  with  a  period  equal  to  the 
period  of  correlation  of  the  process.  In  principle  the  dispersion  may 
be  a  precise  period,  equal  to  any  part  of  the  period  of  correlation, 
and  it  also  may  be  degenerate. 

The  first  assertion  follows  directly  from  formula  (1)  and  the  fact 
that  T  is  the  precise  period  of  covariation.  For  proof  of  the  second 
of  formulas  (3)  and  (4)  we  find 


o*(0=  V 


(5) 


Assume  with  a  fixed  integer  m  only  the  initial  values  of 


KK  (5„p  (0)  0,  p= — oo,  ou)  will  be  different  from  zero.  Then 


from  which,  taking  into  account  that  e®*  =  1 


we  see  that  the 


expression 


L.  is  the  precise  period  of  dispersion,  i.e. 


2 


®  ^  .  Dispersion  is  degenerate  when  Bo(0)j!j(feO,  and  all  the 

remaining  Bj^{0)=0.  These  conditions  are  necessary  and  sufficient  as  a 
result  of  the  uniqueness  of  the  Fourier  expansion. 

Subsequently  we  will  assume  that  PKSP  is  harmonizable  (optionally 
absolutely  in  the  meaning  of  Rozanov  [4])  and  its  two- f r equency 
spectral  functional  is  differentiable  in  the  sense  of  the  theory  of 
generalized  functions,  i.e. 

r(s.i)^J  J  (X.  R)  =  ?  T 

~it  at 

Then  its  derivative  s{  ^  )  is  a  covariation  of  harmonics  with 
frequencies  A.  and  ^  [6]  and  has  the  form 

OO 

s(X.fi)  =  £r(X)Z^=  £  (6) 

Am — go 


where  /k(’)=F'k(')  Z'(*)  -  generalized  derivative  in  the  sense  of 

Ito  [3]  of  the  spectral  measure  of  the  process,  entering  into  its 
representation 


5  (0  -  J  e"*</Z  (X)  -  J 


(7) 


From  the  representation  of  the  function  of  covariation 


(8) 


is  valid  in  the  meaning  of  the  theory  of  generalized  functions  for  any 
harmonizable  SP,  if  [2] 


fiii- V.  v)  =  $(v.  n). 


from  formula  (6)  we  find  that  the  function  of  covariation  of  harmonics 


Figure  1  . 

Definition  3.  We  will  call  the  functions  the  spectral 

components  (SpK).  Then  from  the  last  equality  and  formula  (2)  the 
validity  of  the  theorem  follows  directly. 

Theorem  2.  The  k-th  spectral  component  and  the  k-th  harmonic  of 
change  in  the  spectrum  of  PKSP  coincide. 

Theorem  3.  Each  harmonic  with  a  frequency  Aq  of  expansion  (7)  of 
the  PKSP  is  correlated  only  with  harmonics  of  frequencies 

,  kss—Qo,  oo  ,  and  is  not  correlated  with  the  remainder.  The 

coefficient  of  correlation  of  the  harmonics  and  Aq  +  is 

determined  by  the  value  iAo)  of  the  complex-valued  kg-th  SpK 

f.  {  •)  (Fig.  1)  .  ' 

^  0 

For  establishing  the  bond  between  the  spectral  components  and  the 
positive  and  negative  indexes  we  will  use  the  correlation  [2]  for  KK : 

_ _  «^ii 

Bk  (u)  =  (—  u) «  . 


Then,  taking  the  first  of  formulas  (3)  into  account,  we  find 

The  sets  of  frequencies  "I"  *  form 

partitioning  of  the  axis  of  frequencies  into  classes  of  equivalence, 
since  the  magnitude  of  correlation  depends  not  only  on  k,  but  also  on 

X  0  • 
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Partitioning  will  take  place  only  in  the  case  of  PK  of  white  noise, 
defined  (see  below)  as  a  process  with  covariation  r(s,  t)=g{t)  (5(s-t), 

where  -j  *  ^  -  periodic  function,  since  for  it  in  formula 

(6)  ^k(  A.)cj^,  and  it  possesses  stationarily  correlated  harmonics, 

i.e,  s{  ^  )=s{  A- ^  )  . 


Let  us  note  the  validity  of  the  general  theorem. 


Theorem  4.  We  will  present  the  covariation  of  any  process  which 
is  stationarily  correlated  in  spectrum  in  the  form  of  the  product 
g(t)  S  (s-t),  where  g(t)^0,  and  conversely,  any  nonstationary  noise  is  a 
harmonizable  process  with  stationarily  correlated  harmonics. 

The  last  assertion  follows  directly  from  the  theorem  of  collation 
of  spectra  and  the  Fou r i e r -p resen ta t ion  of  the  5 -function.  Then 


r(s.0^8{nHt-s)^  I 


and  G(*)  -  Four ier- image  of  the  function  g(*)* 

Since  the  function 


•  T 

-jg-j  b(t.  u)dt  ==>  Y  ^ b(t,  u)dt  ^  ^  ^ 

-•  0 

exists  and  is  defined  positively,  then  /q  (  ^  i  describes  the 

powers  of  the  harmoraincs  of  the  components  of  the  process.  On  the 
other  hand,  the  function,  determined  by  the  first  of  equalities  (10), 
is  called  the  covariation  of  a  stationary  approximation  to  a 
nonstationary  process.  Everything  taken  together  shows  that  the 
theorem  is  valid. 

Theorem  5.  A  stationary  approximation  to  PKSP  -  this  is  a 
stationary  process  consisting  of  the  same  harmonic  components  and  with 
the  same  powers  as  PKSP,  but  that  are  not  correlated.  Averaging  (10) 
is  equivalent  to  elimination  of  the  correlation  of  harmonics,  and  the 
power  Bo(0)  of  stationary  approximation  is  equal  to  the  power  of  the 
PKSP  which  is  average  for  the  period. 

Actually,  from  formulas  (10)  and  (4)  it  is  evident  that 

r 

O*(0df  •  In  that  particular  case  when  b{t,  •)  is  almost  a 


5 


periodic  function,  functions  will  be  almost  periodic. 


and  if 


—  OO,  OO) 


-  set  of  indices  of  a  Fourier  function  b(t, 
the  function  of  the  frequency  covariation 


) ,  then 


s(X,  fi)  =» 


o»  e» 


(11) 


is  different  from  zero  only  in  the  points  of  intersection  of  the  set 
|x  =  it —  X,,|  with  the  set  {  ^  .  It  follows  from  here  that  only 


.  Will  be  correlated. 


In 


2.-t 


particular,  for  periodic  processes  [2,6]  ~f 


and  all  the 


harmonics  are  correlated,  and  the  spectrum  is  averaged  in  points 
2.n  .  2.1  ’ 


j  ,  ft.  /  =.  —  OO.  OO 

In  a  general  case  the  following  theorem  takes  place. 


Theorem  6.  If  the  function  b{t,  •)  is  periodic,  then  its  period 
is  a  multiple  of  the  period  of  correlation. 


Assume  the  period  of  the  function  b(t,  •)  is  equal  to  Tq .  Then  on 
the  basis  of  formula  (3)  the  function  possesses  a  Fourier  expansion 

«a  OO  •2;i  f  *  ^  \ 

6(/.  u)  =  V  V  b^,e  ^ ,  from  which  it  is  easy  to  see  that  all  the 

4=»— »  09 

KK  will  be  periodic  with  the  period  Tq: 


— 30 


and  in  this  case  in  conformity  with  formula  (9)  ^  ^ 


X  6 


From  this  expression  and  the  identity  which  is  known 


from  the  theory  of  generalized  functions  0-  S  (OSO  it  follows  that 
when  *^1  ^  t  ft  =  —  coj  =  0  ^  then  f  {  ,  X  )~s0  ,  and 


then  on  the 
formula  (8) 
degenerated 


basis  of  the  uniqueness  of  the  Fourier  transform  from 
we  have  the  identity  b(t,  u )  st  0 ,  i.e.,  the  PKSP  is 

to  zero.  In  the  case  of  non-degeneration  the  inclusion  of 


6 


takes  place,  since  T  is  the  precise  period  of 
correlation.  Since  both  these  sets  of  points  form  arithmetic 

progressions  with  the  differences  -jr  and  ^  respectively,  then 

2.1  2n 

there  exists  such  a  natural  number  p  that  T-=P77  •  The  latter  is 

equivalent  to  the  equality  To=pT,  proving  the  theorem. 

Subsequently  it  will  be  convenient  to  use  the  following 
def initions . 

Definition  4.  In  a  general  representation  of  PKSP  [2] 

EM-  S 

— o® 


where  {?!,(•)}  ~  stationary  stationarily  bound  processes  (stationary 
components  (SK  PKSP). 

De  f i nit  ion  5 .  We  will  designate  as  periodically  correlated  white 
noise  a  PKSP  for  which  in  formula  (6)  •^  (  ^  ) =  ^  foT  ^  ®nd  k. 

From  definition  5  it  is  evident  that 

r(s.  0  =  2n  S  ^  6(s- f)  =  g(0S(s~  0, 


i.e.,  for  PK  of  white  noise  b{t,  u)=g(t)<r  (u),  where  g{‘)^0  -  periodic 
function  of  the  period,  equal  to  the  period  of  correlation. 

Theorem  7.  if  SK  PKSP  -  white  noises,  then  it  degenerates  into  PK 
white  noise  . 

Act!  illy,  using  formula  [2] 


t.i 


-♦(*<)« 


i  ^  iltl-lu) 


(13) 


and  the  condition  of  noncorrelation  of  .^Krj^  j(u)-A]^^j<r(u),  and  also 
the  property  f  (x)  •  (T  (x)=  /  (0)  (  X ) ,  we  obta i n 
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6(/.u)=  S 

ib* — oa 


fi(«)  =  ir(0fi  {«). 


where  g(’)  -  periodic  function  of  period  T.  It  is  non-negative  when 

Al , j>0. 

Theorea  8.  If  SK  are  uncorrelated,  then  PKSP  degenerates  into  a 
stationary  SP,  and  if  they  degenerate  into  random  variables,  then  the 
process  will  be  periodic,  if  these  variables  are  uncorrelated  the 
process  will  be  periodic  stationary  [6]. 

Proof  follows  from  formulas  (9)  and  (10). 


Definition  6.  We  will  designate  as  periodic  white  noise  a  process 


for  which  in  formula  (11)  when 


2n 


Then  its  covariation 


ht  ,  m 


r(3,{) 


V 


which  taking  into  account  the  Poisson  formula  from  the  theory  of 
generalized  functions 


hm  at 


giving  a  representation  of  the  peiodic 


^-function  2 

/-^co 


can  be  written  in  the  form  r(s,  t)=Tg(t)  ^^(t-s),  which  proves  the 
statement . 


Theorem  9.  Periodic  noise  does  not  coincide  with  PK  noise. 


The  following  theorems  are  important  for  an  understanding  of  the 
structure  of  PKSP. 


Theorem  10.  The  state  of  being  correlated  by  stationary  PKSP 
components  and  the  state  of  being  correlated  by  its  harmonic  components 
are  equivalent. 

Since  SK  are  correlated,  and  by  definition  4  they  are  stationarily 
bound,  then  their  mutual  covariation  permits  the  representation 


8 


where  EdZi  (X)  dZj  (|i)  =*  dS/./  (X)d(X  —  /  and  2j^(*)  -  random  measure, 

corresponding  to  the  process  (  ‘  )  in  expansion  (7).  Then  from 
formulas  (3)  and  (13)  we  find 

Urn — C» 

from  which  ^ 

from  which  the  validity  of  the  theorem  follows  directly,  i.e.,  the 
correlated  state  of  the  harmonics  implies  the  correlated  state  of 
stationary  components,  and  vice  versa.  In  particular,  from  here  the 

spectral  function  of  the  statinary  approximation  to  PKSP 

«» 

S/./  —  /  -yj  ,  i.e.,  is  equal  to  the  sum  of  the  displaced 

spectral  functions  of  stationary  components,  and  the  shift  is  equal  to 
the  frequency  of  the  correlated  state  of  PKSP  multiplied  by  the  SK 
number . 

Theorem  11.  An  invariant  ideal  band-pass  filter  separates  from 
the  PKSP  its  uncorrelated  harmonic  components  when,  and  only  when,  the 
width  of  its  bandwidth  is  less  than  the  frequency  of  the  correlated 
state  of  PKSP. 

Proof  follows  from  the  fact  that  the  spectral  density  of  the 
transform^f  the  harmonized  process  with  the  help  of  the  invarxant 
operator  <P  with  the  frequency  characteristic  <t>  {  "X  )  is  expressed 
through  the  spectral  density  of  the  process  by  the  correlation  [2,6] 

s.^(X.fi)  =  (I>(X)'i^Sj{X,ii), 


which  taking  formula  (6)  into  account  takes  the  form 

«•  ■ 

5.,(X.^)  =  a)(X)'  5] 

Assume  now  that  the  frequency  characteristic  of  the  filter  (^)  is 

equal  to  zero  outside  of  the  interval 

i.e.,  <I>x.(>.)-(I>,,(X)xa.(?.)  ,  where  {•)  -  characteristic  function  o'"  the 
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set  A .  Then 


^  fH)  fi  (fi  --  X  +  * 


Considering  the  property  of  the  (J-function 
%^  ( X ) • ^( X ) ^ ( x-a ) a )/{ a ) ,  we  see  that  the  product  of  the  last  three 

M.2x,andasa 

<  y- 

result  of  the  independence  of  A,  and  ^  from  the  inequalities 

|^-Xl<|ji|  +  |Xl  it  turns  out  that  the  product  is  different  from  zero 
only  when  and  |X|<-|S-  ,  Under  these  conditions  the  mentioned 

product  is  equal  to  g 6 (X  —  (i),  ,  and  then 

%  =  XI^  (X)  /„  (X)  I  (X)  P  6  (^  -  X). 

From  here  the  correlation  function  of  the  transform 

(s.  0  ^  1  (D^  (X)  |*dX. 

i.e,,  the  filtration  of  the  PKSP  with  a  bandwidth  less  than  the 
frequency  of  its  correlated  state  is  equivalent  to  the  filtration  of 
its  stationary  approximation.  The  sufficiency  of  the  condition  of  the 
theorem  follows  from  the  fact  that  if  A=[A^-^,  +  A]  is  designated, 

then 

X,4(X)X,,(fi)5(x  =  Xa  {X)6(X-}i)6,, 

only  in  the  case  when  the  set  A  x  A  does  not  intersect  with  one  of  the 

_  . _ .in  .  _ _ 

straight  lines  p  —  A  R  ^  .«=  oo.oo  ^  with  the  exception  of  the 
diagonal  /•  -  A,  (Figure  2). 

Another  proof  of  this  theorem  can  be  obtained  if  in  expansion  (12) 
of  the  PKSP  each  stationary  component  is  replaced  by  its  Cramer  - 
Kolmogorov  expansion 


factors  ditters  trom  zero  only  in  tne  case  wren 


l»{0=  J  <"'‘^*(X). 


Then  the  harmonic  expansion  (7)  for  PKSP  takes  the  form 


t(0=  ]«“<<[  2 

—o*  *«• — am  J 


from  where  follows  directly  the  type  of  correlated  state  of  the 
harmonics  mentioned  in  theorem  3: 

Thus,  in  the  representation  of  PKSP  of  type  (7) 

— cm 

and  in  the  sense  of  the  theorem  of  generalized  functions 


/(v.X)  =  £Z'(v)2'(v-rX). 

Def inition  7 .  We  will  call  a  harmonizable  SP  degenerate  if  its 
dispersion  is  constant. 

According  to  theorem  1  such  a  process  may  be  nonstationary. 

Theorem  12.  A  non-degenerate  harmonizable  SP  with  periodic 
dispersion  belongs  to  the  class  of  almost  periodically  correlated. 

For  proof  of  this  assertion  from  formula  (8)  we  determine  the 
dispersion 


J  f(v,  X)dv(iX. 


(14) 


it  is  evident  from  expression  (14)  that  dispersion  can  be  periodic 
with  a  precise  period  Tg  only  when 


The  latter  in  turn  is  possible  only  when 


/(V.  X)  =  /j(v.  X)  +  /,(v.  X)  £  6(X-X,)  =  /,(v.  X)-l-  V  /^(v)6(X-X,). 


where  7^{V)  =  )y{  A  )•  Iri  this  case  the  following  conditions  should 
be  fulfilled. 


From  the  last  equality  and  formula  (14)  it  follows  that  dispersion 
of  the  component  of  a  random  process,  possessing  a  continuous  spectrum 
A  )/  is  equal  to  zero.  Since  the  process  is  centered,  then 
this  component  is  identically  equal  to  zero,  which  also  proves  the 
theorem . 

It  follows  from  expression  (14)  that  if  z  is  taken  as  yq,  then  the 
inequality  (12)  will  be  fulfilled.  For  a  more  rapid  convergence  of 
formula  (11)  when  y>l  it  is  possible  to  use  as  the  initial 
approximation  the  magnitude 
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The  table  gives  data  on  the  absolute  error  of  formulas  (6)  and 
(11),  and  also  of  the  iteration  formula  of  the  third  order 


y,+,  =  y, — ^  (y,)  ~  *1  +  t  y<  I*!*  (y)  - 


.*  a 


(16) 


in  the  case  of  the  initial  approximation  (13)  and 


yo=*j  = 


z  when  0  <  x  <  0.82; 

,  ^  r  >  0.82. 

(  z 1  when 


(17) 
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